Adopting duplex sequencing technology for genetic toxicity testing: A proof-of-concept mutagenesis experiment with N-ethyl-N-nitrosourea (ENU)-exposed rats.

Duplex sequencing (DS) is an error-corrected next-generation sequencing method in which molecular barcodes informatically link PCR-copies back to their source DNA strands, enabling computational removal of errors in consensus sequences. The resulting background of less than one artifactual mutation per 107 nucleotides allows for direct detection of somatic mutations. TwinStrand Biosciences, Inc. has developed a DS-based mutagenesis assay to sample the rat genome, which can be applied to genetic toxicity testing. To evaluate this assay for early detection of mutagenesis, a time-course study was conducted using male Hsd:Sprague Dawley SD rats (3 per group) administered a single dose of 40 mg/kg N-ethyl-N-nitrosourea (ENU) via gavage, with mutation frequency (MF) and spectrum analyzed in stomach, bone marrow, blood, and liver tissues at 3 h, 24 h, 7 d, and 28 d post-exposure. Significant increases in MF were observed in ENU-exposed rats as early as 24 h for stomach (site of contact) and bone marrow (a highly proliferative tissue) and at 7 d for liver and blood. The canonical, mutational signature of ENU was established by 7 d post-exposure in all four tissues. Interlaboratory analysis of a subset of samples from different tissues and time points demonstrated remarkable reproducibility for both MF and spectrum. These results demonstrate that MF and spectrum can be evaluated successfully by directly sequencing targeted regions of DNA obtained from various tissues⁠, a considerable advancement compared to currently used in vivo gene mutation assays.

Adopting Duplex Sequencing™ Technology for Genetic Toxicity Testing: A Proof-of-Concept Mutagenesis Experiment with N-Ethyl-N-Nitrosourea (ENU)-Exposed Rats.

Duplex sequencing (DuplexSeq) is an error-corrected next-generation sequencing (ecNGS) method in which molecular barcodes informatically link PCR-copies back to their source DNA strands, enabling computational removal of errors by comparing grouped strand sequencing reads. The resulting background of less than one artifactual mutation per 10 7 nucleotides allows for direct detection of somatic mutations. TwinStrand Biosciences, Inc. has developed a DuplexSeq-based mutagenesis assay to sample the rat genome, which can be applied to genetic toxicity testing. To evaluate this assay for early detection of mutagenesis, a time-course study was conducted using male Hsd:Sprague Dawley SD rats (3 per group) administered a single dose of 40 mg/kg N-ethyl-N-nitrosourea (ENU) via gavage, with mutation frequency (MF) and spectrum analyzed in stomach, bone marrow, blood, and liver tissues at 3 h, 24 h, 7 d, and 28 d post-exposure. Significant increases in MF were observed in ENU-exposed rats as early as 24 h for stomach (site of contact) and bone marrow (a highly proliferative tissue) and at 7 d for liver and blood. The canonical, mutational signature of ENU was established by 7 d post-exposure in all four tissues. Interlaboratory analysis of a subset of samples from different tissues and time points demonstrated remarkable reproducibility for both MF and spectrum. These results demonstrate that MF and spectrum can be evaluated successfully by directly sequencing targeted regions of DNA obtained from various tissues, a considerable advancement compared to currently used in vivo gene mutation assays. HIGHLIGHTS: DuplexSeq is an ultra-accurate NGS technology that directly quantifies mutationsENU-dependent mutagenesis was detected 24 h post-exposure in proliferative tissuesMultiple tissues exhibited the canonical ENU mutation spectrum 7 d after exposureResults obtained with DuplexSeq were highly concordant between laboratoriesThe Rat-50 Mutagenesis Assay is promising for applications in genetic toxicology.

Genetic toxicity testing using human in vitro organotypic airway cultures: Assessing DNA damage with the CometChip and mutagenesis by Duplex Sequencing.

The organotypic human air-liquid-interface (ALI) airway tissue model has been used as an in vitro cell culture system for evaluating the toxicity of inhaled substances. ALI airway cultures are highly differentiated, which has made it challenging to evaluate genetic toxicology endpoints. In the current study, we assayed DNA damage with the high-throughput CometChip assay and quantified mutagenesis with Duplex Sequencing, an error-corrected next-generation sequencing method capable of detecting a single mutation per 107 base pairs. Fully differentiated human ALI airway cultures were treated from the basolateral side with 6.25 to 100 µg/mL ethyl methanesulfonate (EMS) over a period of 28 days. CometChip assays were conducted after 3 and 28 days of treatment, and Duplex Sequencing after 28 days of treatment. Treating the airway cultures with EMS resulted in time- and concentration-dependent increases in DNA damage and a concentration-dependent increase in mutant frequency. The mutations observed in the EMS-treated cultures were predominantly C → T transitions and exhibited a unique trinucleotide signature relative to the negative control. Measurement of physiological endpoints indicated that the EMS treatments had no effect on anti-p63-positive basal cell frequency, but produced concentration-responsive increases in cytotoxicity and perturbations in cell morphology, along with concentration-responsive decreases in culture viability, goblet cell and anti-Ki67-positive proliferating cell frequency, cilia beating frequency, and mucin secretion. The results indicate that a unified 28-day study can be used to measure several important safety endpoints in physiologically relevant human in vitro ALI airway cultures, including DNA damage, mutagenicity, and tissue-specific general toxicity.

Direct quantification of in vivo mutagenesis and carcinogenesis using duplex sequencing.

The ability to accurately measure mutations is critical for basic research and identifying potential drug and chemical carcinogens. Current methods for in vivo quantification of mutagenesis are limited because they rely on transgenic rodent systems that are low-throughput, expensive, prolonged, and do not fully represent other species such as humans. Next-generation sequencing (NGS) is a conceptually attractive alternative for detecting mutations in the DNA of any organism; however, the limit of resolution for standard NGS is poor. Technical error rates (∼1 × 10-3) of NGS obscure the true abundance of somatic mutations, which can exist at per-nucleotide frequencies ≤1 × 10-7 Using duplex sequencing, an extremely accurate error-corrected NGS (ecNGS) technology, we were able to detect mutations induced by three carcinogens in five tissues of two strains of mice within 31 d following exposure. We observed a strong correlation between mutation induction measured by duplex sequencing and the gold-standard transgenic rodent mutation assay. We identified exposure-specific mutation spectra of each compound through trinucleotide patterns of base substitution. We observed variation in mutation susceptibility by genomic region, as well as by DNA strand. We also identified a primordial marker of carcinogenesis in a cancer-predisposed strain of mice, as evidenced by clonal expansions of cells carrying an activated oncogene, less than a month after carcinogen exposure. These findings demonstrate that ecNGS is a powerful method for sensitively detecting and characterizing mutagenesis and the early clonal evolutionary hallmarks of carcinogenesis. Duplex sequencing can be broadly applied to basic mutational research, regulatory safety testing, and emerging clinical applications.

Ultra-accurate Duplex Sequencing for the assessment of pretreatment ABL1 kinase domain mutations in Ph+ ALL.

Mutations of ABL1 are the dominant mechanism of relapse in Philadelphia chromosome-positive acute lymphoblastic leukemia (Ph + ALL). We performed highly accurate Duplex Sequencing of exons 4-10 of ABL1 on bone marrow or peripheral blood samples from 63 adult patients with previously untreated Ph + ALL who received induction with intensive chemotherapy plus a BCR-ABL1 TKI. We identified ABL1 mutations prior to BCR-ABL1 TKI exposure in 78% of patients. However, these mutations were generally present at extremely low levels (median variant allelic frequency 0.008% [range, 0.004%-3.71%] and did not clonally expand and lead to relapse in any patient, even when the pretreatment mutation was known to confer resistance to the TKI received. In relapse samples harboring a TKI-resistant ABL1 mutation, the corresponding mutation could not be detected pretreatment, despite validated sequencing sensitivity of Duplex Sequencing down to 0.005%. In samples under the selective pressure of ongoing TKI therapy, we detected low-level, emerging resistance mutations up to 5 months prior to relapse. These findings suggest that pretreatment ABL1 mutation assessment should not guide upfront TKI selection in Ph + ALL, although serial testing while on TKI therapy may allow for early detection of clinically actionable resistant clones.

Ultra-Sensitive TP53 Sequencing for Cancer Detection Reveals Progressive Clonal Selection in Normal Tissue over a Century of Human Lifespan.

High-accuracy next-generation DNA sequencing promises a paradigm shift in early cancer detection by enabling the identification of mutant cancer molecules in minimally invasive body fluid samples. We demonstrate 80% sensitivity for ovarian cancer detection using ultra-accurate Duplex Sequencing to identify TP53 mutations in uterine lavage. However, in addition to tumor DNA, we also detect low-frequency TP53 mutations in nearly all lavages from women with and without cancer. These mutations increase with age and share the selection traits of clonal TP53 mutations commonly found in human tumors. We show that low-frequency TP53 mutations exist in multiple healthy tissues, from newborn to centenarian, and progressively increase in abundance and pathogenicity with older age across tissue types. Our results illustrate that subclonal cancer evolutionary processes are a ubiquitous part of normal human aging, and great care must be taken to distinguish tumor-derived from age-associated mutations in high-sensitivity clinical cancer diagnostics.

dNTP pool levels modulate mutator phenotypes of error-prone DNA polymerase ε variants.

Mutator phenotypes create genetic diversity that fuels tumor evolution. DNA polymerase (Pol) ε mediates leading strand DNA replication. Proofreading defects in this enzyme drive a number of human malignancies. Here, using budding yeast, we show that mutator variants of Pol ε depend on damage uninducible (Dun)1, an S-phase checkpoint kinase that maintains dNTP levels during a normal cell cycle and up-regulates dNTP synthesis upon checkpoint activation. Deletion of DUN1 (dun1Δ) suppresses the mutator phenotype of pol2-4 (encoding Pol ε proofreading deficiency) and is synthetically lethal with pol2-M644G (encoding altered Pol ε base selectivity). Although pol2-4 cells cycle normally, pol2-M644G cells progress slowly through S-phase. The pol2-M644G cells tolerate deletions of mediator of the replication checkpoint (MRC) 1 (mrc1Δ) and radiation sensitive (Rad) 9 (rad9Δ), which encode mediators of checkpoint responses to replication stress and DNA damage, respectively. The pol2-M644G mutator phenotype is partially suppressed by mrc1Δ but not rad9Δ; neither deletion suppresses the pol2-4 mutator phenotype. Thus, checkpoint activation augments the Dun1 effect on replication fidelity but is not required for it. Deletions of genes encoding key Dun1 targets that negatively regulate dNTP synthesis, suppress the dun1Δ pol2-M644G synthetic lethality and restore the mutator phenotype of pol2-4 in dun1Δ cells. DUN1 pol2-M644G cells have constitutively high dNTP levels, consistent with checkpoint activation. In contrast, pol2-4 and POL2 cells have similar dNTP levels, which decline in the absence of Dun1 and rise in the absence of the negative regulators of dNTP synthesis. Thus, dNTP pool levels correlate with Pol ε mutator severity, suggesting that treatments targeting dNTP pools could modulate mutator phenotypes for therapy.

Emergence of DNA polymerase ε antimutators that escape error-induced extinction in yeast.

DNA polymerases (Pols) ε and δ perform the bulk of yeast leading- and lagging-strand DNA synthesis. Both Pols possess intrinsic proofreading exonucleases that edit errors during polymerization. Rare errors that elude proofreading are extended into duplex DNA and excised by the mismatch repair (MMR) system. Strains that lack Pol proofreading or MMR exhibit a 10- to 100-fold increase in spontaneous mutation rate (mutator phenotype), and inactivation of both Pol δ proofreading (pol3-01) and MMR is lethal due to replication error-induced extinction (EEX). It is unclear whether a similar synthetic lethal relationship exists between defects in Pol ε proofreading (pol2-4) and MMR. Using a plasmid-shuffling strategy in haploid Saccharomyces cerevisiae, we observed synthetic lethality of pol2-4 with alleles that completely abrogate MMR (msh2Δ, mlh1Δ, msh3Δ msh6Δ, or pms1Δ mlh3Δ) but not with partial MMR loss (msh3Δ, msh6Δ, pms1Δ, or mlh3Δ), indicating that high levels of unrepaired Pol ε errors drive extinction. However, variants that escape this error-induced extinction (eex mutants) frequently emerged. Five percent of pol2-4 msh2Δ eex mutants encoded second-site changes in Pol ε that reduced the pol2-4 mutator phenotype between 3- and 23-fold. The remaining eex alleles were extragenic to pol2-4. The locations of antimutator amino-acid changes in Pol ε and their effects on mutation spectra suggest multiple mechanisms of mutator suppression. Our data indicate that unrepaired leading- and lagging-strand polymerase errors drive extinction within a few cell divisions and suggest that there are polymerase-specific pathways of mutator suppression. The prevalence of suppressors extragenic to the Pol ε gene suggests that factors in addition to proofreading and MMR influence leading-strand DNA replication fidelity.

Mutator suppression and escape from replication error-induced extinction in yeast.

Cells rely on a network of conserved pathways to govern DNA replication fidelity. Loss of polymerase proofreading or mismatch repair elevates spontaneous mutation and facilitates cellular adaptation. However, double mutants are inviable, suggesting that extreme mutation rates exceed an error threshold. Here we combine alleles that affect DNA polymerase δ (Pol δ) proofreading and mismatch repair to define the maximal error rate in haploid yeast and to characterize genetic suppressors of mutator phenotypes. We show that populations tolerate mutation rates 1,000-fold above wild-type levels but collapse when the rate exceeds 10⁻³ inactivating mutations per gene per cell division. Variants that escape this error-induced extinction (eex) rapidly emerge from mutator clones. One-third of the escape mutants result from second-site changes in Pol δ that suppress the proofreading-deficient phenotype, while two-thirds are extragenic. The structural locations of the Pol δ changes suggest multiple antimutator mechanisms. Our studies reveal the transient nature of eukaryotic mutators and show that mutator phenotypes are readily suppressed by genetic adaptation. This has implications for the role of mutator phenotypes in cancer.

Antimutator variants of DNA polymerases.

Evolution balances DNA replication speed and accuracy to optimize replicative fitness and genetic stability. There is no selective pressure to improve DNA replication fidelity beyond the background mutation rate from other sources, such as DNA damage. However, DNA polymerases remain amenable to amino acid substitutions that lower intrinsic error rates. Here, we review these 'antimutagenic' changes in DNA polymerases and discuss what they reveal about mechanisms of replication fidelity. Pioneering studies with bacteriophage T4 DNA polymerase (T4 Pol) established the paradigm that antimutator amino acid substitutions reduce replication errors by increasing proofreading efficiency at the expense of polymerase processivity. The discoveries of antimutator substitutions in proofreading-deficient 'mutator' derivatives of bacterial Pols I and III and yeast Pol δ suggest there must be additional antimutagenic mechanisms. Remarkably, many of the affected amino acid positions from Pol I, Pol III, and Pol δ are similar to the original T4 Pol substitutions. The locations of antimutator substitutions within DNA polymerase structures suggest that they may increase nucleotide selectivity and/or promote dissociation of primer termini from polymerases poised for misincorporation, leading to expulsion of incorrect nucleotides. If misincorporation occurs, enhanced primer dissociation from polymerase domains may improve proofreading in cis by an intrinsic exonuclease or in trans by alternate cellular proofreading activities. Together, these studies reveal that natural selection can readily restore replication error rates to sustainable levels following an adaptive mutator phenotype.

Glucosidase-beta variations and Lewy body disorders.

It has been proposed that there is an increased frequency of glucosidase-beta mutations in Lewy body disorders. Our comprehensive DNA sequencing approach found a small number of glucosidase-beta mutations in 101 neuropathologically defined Lewy body disease cases (3%) compared to 99 healthy post-mortem controls (1%); odds ratio 3.0 (95% CI: 0.3-29, p=0.3). All three affected carriers were classified as diffuse Lewy body disease (n=3/50; 6%). Our study suggests glucosidase-beta variants have a limited role in susceptibility to Lewy body disease in North America.

Linkage disequilibrium and association analysis of alpha-synuclein and alcohol and drug dependence in two American Indian populations.

BACKGROUND: Alpha-synuclein is involved in dopaminergic neurotransmission and has been implicated in a number of neurodegenerative disorders, such as Parkinson's disease. Recent studies, in humans and in rat and monkey models, have suggested that alpha-synuclein may play a role in the development and maintenance of certain addictive disorders. METHODS: Fifteen single-nucleotide polymorphisms (SNPs) in the alpha-synuclein gene (SNCA) and 1 upstream microsatellite repeat (NACP-REP1) were assayed in Southwest (SW; n=514) and Plains (n=420) American Indian populations. Patterns of linkage disequilibrium (LD) at SNCA were determined for the 2 populations and compared with Caucasian, African, and Asian populations in the HapMap database (http://www.hapmap.org). Assayed alleles and constructed haplotypes in the study populations were tested for association with 4 clinical phenotypes [alcohol dependence, alcohol use disorders, drug dependence, and drug use disorders (lifetime diagnoses)] as well as with 2 symptom count phenotypes (all 18 questions and the 8 questions diagnostic for alcohol dependence). RESULTS: Patterns of LD at SNCA were similar in both Indian populations and were consistent with the LD structure in other populations as reflected in the HapMap database. Single allele tests revealed significant associations between 4 SNPs and drug dependence in the SW population and between 2 of those SNPs plus 2 other SNPs and drug dependence in SW males only. In the Plains population, a significant association was detected only in males between 2 SNPs and alcohol use disorders and between 1 SNP and alcohol dependence. In the SW population, 1 SNP was marginally significant with the total symptom count. However, in all cases, the support was modest and disappeared with correction for multiple comparisons. No association was found between constructed haplotypes and any of the phenotypes in either population. CONCLUSIONS: Despite modest support for association between multiple SNCA SNPs and several of the addictive disorders tested in this study, statistical significance disappeared after correction for multiple testing. Thus, our data do not support a role for a variant in the SNCA gene that contributes to alcohol or drug addiction in the 2 studied American Indian populations. Future research may focus on variants in the promoter region that could cause the changes in mRNA and protein levels observed in previous studies.